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Abstract 

An approach is developed, which essentially simplifies derivation of kinetic equation for polar- 
ization in non-relativistic scattering. This approach is applicable for collision of projectile particles 
with a target for any spins of colliding particles. The most detailed consideration is performed for 
the case of spin 1/2 projectile particle interacting with spin 1/2 particle of a target. The solution 
of the kinetic equation for the case of zero initial polarization is presented. 
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I. INTRODUCTION 



It is known that the collision of high energy beams of polarized particles gives additional 
very important information as compared with the case of unpolarized particles. However, 
;o obtain high energy beam of heavy particles with noticeable polarization is hard problem 
l sj. Besides, theoretical description of the kinetics of polarization is also nontrivial even 
in non-relativistic case and needs its further development, see Review (4]. In Ref. J5J, the 
kinetics of the polarization buildup during the interaction of stored protons (antiprotons) 
with a polarized target was investigated. The kinetic equation was written in terms of 
;he cross section corresponding to spin flip transition and spin non-flip transition. In Ref. 
5], the quantization axes was directed along the polarization vector £<r of the target. If 
£t is parallel or perpendicular to the momentum of particles in the beam, and the initial 
polarization of these particles is zero, then it follows from the arguments of parity that 
the polarization, arising as a result of interaction, is directed along Similar approach 
is used in Ref. (6| at the consideration of polarization effects in non-relativistic electron- 
proton scattering. If is n °t parallel or perpendicular to the momentum of particles in 
the beam, then it is necessary to take into account additional terms in the kinetic equation 
which lead to rotation of polarization vector in the process of the polarization buildup 
The kinetic equation for the density matrix was derived for the first time in Ref. (8j by 
solving the quantum Liouville equation and expressing its solution via a spin-dependent 
scattering amplitude. The equation obtained in Ref. |8[ was applied to the spin-exchange 
optical pumping, and it was shown that the spin-exchange collisions lead to rotation of 
electron spin. The same kinetic equation was later re-derived in Ref. [9| and applied to the 



analysis of the spin evolution, see also 



In the present paper we develop an approach which essentially simplifies derivation of 
kinetic equation describing spin evolution for any spin of the target and particles in the 
beam. We consider the non-relativistic scattering and restrict ourselves for simplicity to 
the case of projectiles remaining in the beam after scattering. This case corresponds to 
small-angle scattering of the projectile. For instance, our approach can be directly applied 
to interaction of heavy particle with the flow of polarized electrons. However, our approach 
can be easily generalized to more complicated cases. 
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II. KINETIC EQUATION 



Let us consider a particle with the spin Si interacting with a flow of polarized particles 
with the spin S^. Below this flow is referred to as a beam. We are going to investigate time 
evolution of < Si >. We assume for simplicity that the mass Mi of the first particle is much 
larger than the mass M2 of the second particle. In this case, in the rest frame of the first 
particle it is possible to consider the first particle as a source of some potential depending on 
the spin operators Si and S^. The operator O is some operator constructed from the spin 
operators, so that [O, r] = and [0,p] = 0, where r is the relative coordinate vector and 
p is the corresponding momentum. The wave function, which corresponds to the scattering 
problem, has the asymptotic form at large distances, 

n ikr 



e lkr + F 

r 



X1X2, (1) 



where k is the initial momentum, xi an d X2 are the spin wave functions of the corresponding 
particles, the operator F depends on no = k/k, n = r/r, and the spin operators. To describe 
kinetics of polarization, we should use the wave packet normalized to unity. We introduce 
the wave function 

^ k (r) = -^=e- Xr Mr), (2) 



where V is some normalization volume, and A is some real parameter which we tend to zero 
at the end of calculations. The factor e~ Xr allows one to perform integration by parts in 
the matrix elements (to use hermiticity of the operators), the system of units h = 1 is used. 
Then we have the usual equation 

j t J drm + k {r)0 H ^ k {r)=iN J dr^+(r)[H,0 H }^ k (r), (3) 

where Oh = e tHt Oe~ lHt is the Heisenberg operator, H is the Hamiltonian of the system, 
iV = 1/V is the density, and [a, b] stands for the commutator of the operators a and b. Using 
the relation Hipk(r) = Eipk(f), we can write Eq. (j3J) as follows, 

j t J drn(r)0 H * k (r) 

= iN J dri>+{r){[e- Xr ,H]0 H e- Xr + e- Xr O H [e- Xr ,H}}Mr)- (4) 

The commutator [e~ Xr , H] is proportional to the small parameter A, and only the contribu- 
tion of large distances r ~ 1/A in the integral over r can compensate this small parameter. 



Therefore, at calculation of the commutator we can leave only the kinetic energy operator 
in the Hamiltonian, 



-Ar 



-Ar V 



iX 



(p 



ne 



-Ar 



e- Xr n ■ p) 



(5) 



' 2M J 2M 

where M « M 2 is the reduced mass. We can also use the asymptotic form (CQ) of the wave 
function. Besides, in the term corresponding to interference of the plane wave and the 
spherical wave, the main contribution to the matrix element is given by the small angle 9 
between vectors n and k, namely 9 2 ~ 1/kr ~ X/k. Finally we obtain the kinetic equation 
d 



dt 



< O >= vNSp < pit) 



2ixi 

dn n F + OF - — (F + (0)O - OF(0) 



(6) 



Here v = k/M, dfl n is the differential of the solid angle corresponding to vector n, F(0) 
is the operator F calculated at n = n , and p(t) is the density matrix which describes the 
spin state of the system, trace is taken over spin indexes of both particles. The density 
matrix equals to pit) = px(t)p2, where pi(t) is the time dependent density matrix of the 
first particle and P2 is the time independent density matrix of the beam. The matrix pi(t) 
should be found as a result of solution of kinetic equation (see below). In (jBJ) we assume 
that .F(O) is finite quantity, otherwise it is necessary to introduce a regularization. 
If we set = 1 then we obtain the unitarity relation 

2m 



Sp p(t) 



dn n F + F-— (F+(0)-F(0) 







(7) 



Eqs. ()6]) and (j7]) are valid for arbitrary spins S% and S 2 of the particles. 



III. THE CASE OF A PARTICLE WITH Si = 1/2. 



Let St = 1/2. Then it follows from Eq. p for O = cr x that 



j t Ci = vNSp{pit) 



2,71 j / 

dQnF+tnF - — (f + (o) CTi - <nF(o) 



Pl(t) = -[l + &(t)-tT 1 ] , Ci(t)=<^> 



(8) 



where <j\ are the Pauli matrices acting on the spin variables of the first particle. The 
unitarity relation (J7J) gives 

~2ni 



Sp 
Sp 



P2 / dn n F + F 



Sp 



P2 (t 1 / dn n F + F 



Sp 



^p 2 (V(0)-F(0)) 
— P2(Tl [F + (0)-FiO)) 



(9) 



For a particle with the initial polarization ( in a moment t, and the polarization 
measured by a detector, the cross section a has the form 



o = Sp 



P (t) j dn n F + Pf F 



Pf = 2 t 1 + ' °il ' 



A = Sp 



/ dQ n F + F 



p(t) / dn n F + (T X F 



(10) 



As a result of scattering the polarization becomes equal to £' = B/A. Performing summation 
over £f and using the unitarity relation (J7J), we write the total cross section a to t as 



Otat = A = A + B ■ Ci , 

A) = ^Sp[p 2 (V(0)-F(0) 

In terms of B and B , Eq. ((HI) reads, 



-Bo = ^-Sp 



P2ffi(F + (0)-F(0) 



(11) 



p 2 (Ci-o-i)(i ?+ (0)o- 1 - 0-^(0) 



(12) 



It is convenient to write Eq. (jH|) also in another form, 

: Ci=vNSp{^p(t) J dQ^F+^^F + F+^F] 
+ ^P 2 [CiX<t 1 ](f+(0) + F(0))}. 

The quantity F can be written as 

F = F + a i ■ F\ , 



(13) 



(14) 



where Fo and Fi are operators acting on the spin variables of the second particle. Then we 
obtain, 



d 
~dt 



Ci = vNSp 2 lp 2 / dQ n F+(F 1 ■ Ci) + (F+ ■ Ci)F 1 - 2(F+ ■ F X )G 



- F+F ) x Ci] - 2z[F+ x Fj] - ^ 2 [(F+(0) + F^O)) x &]} , (15) 

where Sp 2 stands for trace over spin variables of particles from the beam (flow of polarized 
particles). 



5 



A. The case S 2 = 0. 



For So = we have, 



F = fo + Zio-i ■ v 



fo, F x = f x u, u=[nxn ] 



(16) 



where /o and fx are some functions of x = n ■ n . For S 2 = 0, the unitarity relations 
reduces to one nontrivial relation 

f 4-7T 

J dn n [\f \ 2 + u 2 \f 1 \ 2 ] = T lmf (0). (17) 
Using Eq. (fl5|) we arrive at the following equation describing spin relaxation, 

— Ci = -w[Ci + (Ci ■ n )n Q ] , 

u = vN J dn n v 2 \f x \ 2 . (18) 

The solution of this equation reads 

Ci{t) = [Ci(0) - n (Ci(0) • n )]e^ + n (Ci(0) • n ) e - 2 ^ . (19) 

Thus, during relaxation we have not only diminishing of G(£) but also rotation of the 
direction of Ci(t). 



B. The case S 2 = 1/2 . 



For S 2 = 1/2, 



P2 



[1 + C2 • o- 2 ] , f = f + (fxcrx + f 2 a 2 ) ■ v + rvy 2 , 



+ ^ = /i^ + TV 2 . (20) 

Here cr 2 are the Pauli matrices acting on the spin variables of the particle from the beam, £2 
is the time independent polarization of the beam. The functions, fo, fx, and f 2 depend on 
x = n ■ n, and the symmetric tensor T* J is constructed from the vectors n and n. Using 
these definitions, we obtain from Eq. f JT5|) . 

-Q = RV{{ + [Cx X + 

2vN J dtt n |/ X | 2 (z/V - v 2 5 ij ) + Re(T ia *T ja ) - T ab *T ab 5 tJ 



d_ 
~dt 

R ij 



T 



g i = 2vN J dtt n e ijk e abc T ja *T kb Q . 



h _ Re T-(0)C 2 a } , 



(21) 



Here we use the relations, 



dVl n X (n ■ n )T^u k = , / d£l n e abc \m(T ia *T 3b ) = , / dQ n lm(T ta *T tb ) = 



valid for any function x( x )- These relations can be easily proved using the representation 



of the tensor T u 



12| 



T ij = 5 ij f 3 + (nV + nini)f 4 + (nW + n*n J ')/ 5 , 



(22) 



where 73,4,5 are some functions of n ■ n . For 5* 2 = 1/2, the unitarity relations (j5J) reduce to 
two nontrivial relations 



J dQ n [\f \ 2 + z/ 2 |/i| 2 + ^l/ 2 | 2 + T ab *T a& ] = ylm/o(0) 
dtt n {2Re(f*T ab ) + 2Re{flf 2 )v a v b - T ij *T a ?e iaa e jt 
Then we obtain the form of the tensor and the vectors J- and Q 



^ImT aft (0) . (23) 
k 



T = A2C2 + B 2 (C2 ■ n )n , 
g = A 3 C 2 + B 3 {£ 2 ■ n )n , 

where A$ and Bi are some numbers. These numbers are expressed as 

|/t| 2 i/ 2 + r a *r a + niniT ia *T ja 
IM V + T ia *T ia - 3niniT ia *T ja 

/ *fT--3nX T^+/ 2 7i^ 2 



(24) 



A 1 = -vN J da 
B 1 = -vN J dQ 
A 2 = vN^ J dVt n \m 
B 2 = -vN^ J dtt n Im 
A 3 = -2vN [ da 



2tt 

+ x 

2tt, 



r(o)-nX^(o) 



+ —Re 

k 



T"(0) -3nX T ij '(0) 



n\nlT a *T 3a - Re fi^jPT" 



3n l n J T m *T^ a - 3Re\^} ) n 3 T t3 T aa * j + T n T n * - T 3 T 

Note that A x < and A + fi x < 0. 

If £2 = then J 7 = Q = 0, and we obtain from Eq. (J2TJ) 



-12J* 



(25) 



Ci(t) = [Ci(0) - n (Ci(0) ■ n )]e^ + n (Ci(0) • f»o)e<* +BO * , 



(26) 



so that we have depolarization. The component of £1 parallel to no and the component of 
Ci transverse to this vector diminish with different rates. 

Let C2 ^ but [£ 2 x n ] = (£2 is parallel to n ). In this case it is easy to find that 

&(*) = e Mt {cos(u;t)Ci(0) + sin(wt)[Ci(0) x n ]} 
+ [e (Al+Bl) * - e Mi cos(u;t)] (Ci(0) • n ) n 



+ [e(^-l]^±f C2, uj = (A 2 + B 2 )C 2 
For (£ 2 • ^0) = (£2 is perpendicular to n ), the solution reads 

,(Ai+Bi/2)t cos (Q t )^(0) 



(27) 



Ci(«) 



3 Ait _ (Ai+Bi/2)t 



COS 



(m)] 



(Ci(0)-C 2 ) 



Aa 



c 2 



+ - 1] ^ 



+e 



4 



(Ai+Bi/2)t 



sin(f2i) 



A 2 [Ci(0) x C2] 



/> \ (Ci(°) • [Ca x wo]) r > 
(Ci(0) • n ) n - 1 2 - [C2 x n ] 

C2 



(28) 



For the general case, [£2 x n ] 7^ and (£ 2 ■ n ) 7^ 0, it is convenient to write Ci{t) as 



Ci(t) = (a + 7)(Ca • n )n - /3(C 2 • n )[C 2 x n ] - 7C2 , 
where a, (3, and 7 are some functions of time. For we have 

C 2 (t) = « 2 (C 2 • n ) 2 + /3 2 (C 2 • n ) 2 [C 2 x n ] 2 + 7% X n ] 2 



(29) 



(30) 



From Eqs. (J2ID and (J2U) we find 

da 

d/3 



(Ax + B x )a - A 2 [C 2 x n ] 2 /3 + A 3 + B 3 



-£ = A 2 a + A^ + (A 2 + S 2 )7 



c?7 
dt 



-{A 2 + B 2 )(C 2 ■n ) 2 f3 + A ll - A 3 



(31) 



Let us write these equations in the matrix form, 









^ 3 + ^ 













(32) 




V 7 / 




{ ~A 3 J 





The formal solution of this equations has the form 

m = e Ut [m + U-^]-U-^. (33) 

The eigenvalues A» of the matrix U read A, = A\ — B\\i, where 

Af + A 4 2 + (c 1 + c 2 )A, + c 2 = 0, Cl = ^|[C 2 xn ] 2 , c 2 = ^±^l!(C 2 ■ n ) 2 . (34) 

It is possible to show that all A* have negative real parts so that the asymptotic form of 
ip(t) at large t is ip(t — > oo) = —U^ 1 ^, or 

a = + B 2 )[(A 2 + B 2 )(A 3 + B 3 )(C 2 • n ) 2 + A 2 A 3 [C 2 x n ] 2 ] + A 2 (A 3 + 5 3 )} , 

/3 = ^[Ai^aSg - A 3 B 2 ) - B 1 A 3 (A 2 + B 2 )} , 

7 = ^[^2(A 2 + B 2 )(A 3 + B 3 )(C2 ■ n ) 2 + A 3 A 2 2 [C 2 x n ] 2 + A 1 A 3 (A 1 + fix)] , 
W={Ai + B X ){A 2 + 5 2 ) 2 (C 2 ■ n ) 2 + A^fo x n ] 2 + A 2 (Ai + B x ) . (35) 

As should be, this asymptotic form is independent of the initial condition. It is also time 
independent though Ci{t) changes its direction in a process of polarization at finite t. 

IV. CONCLUSIONS 

We have developed an approach which essentially simplifies derivation of the kinetic 
equation describing spin evolution for any spin of the target and particles in the beam. As 
an example, we have considered the non-relativistic scattering of the flow of polarized light 
particles on the heavy particle. For Si = S 2 = 1/2, we have obtained the explicit solution of 
the kinetic equation (Bill , which is expressed via a few constants, Eq. (|25l) . The asymptotic 
form of the solution is independent of the initial condition and time (there is no rotation 
at large time) though Ci(i) changes its direction in a process of polarization. Our approach 
can be easily generalized to more complicated cases. 
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